// This stacks enahos so we have the right survey weights in 2019 quarters and 2020 quarters

cap program drop addSectorLabels
program define addSectorLabels
	** Sectores Completos a CIIU Rev.4
	* Actividad económica
	gen sector1r4=.
	replace sector1r4=1 if p506r4>=100 & p506r4<300
	replace sector1r4=2 if p506r4>=300 & p506r4<500
	replace sector1r4=3 if p506r4>=500 & p506r4<1000
	replace sector1r4=4 if p506r4>=1000 & p506r4<3500
	replace sector1r4=5 if p506r4>=3500 & p506r4<3600
	replace sector1r4=6 if p506r4>=3600 & p506r4<4100
	replace sector1r4=7 if p506r4>=4100 & p506r4<4500
	replace sector1r4=8 if p506r4>=4500 & p506r4<4900
	replace sector1r4=9 if p506r4>=4900 & p506r4<5500
	replace sector1r4=10 if p506r4>=5500 & p506r4<5800
	replace sector1r4=11 if p506r4>=5800 & p506r4<6400
	replace sector1r4=12 if p506r4>=6400 & p506r4<6800
	replace sector1r4=13 if p506r4>=6800 & p506r4<6900
	replace sector1r4=14 if p506r4>=6900 & p506r4<7700
	replace sector1r4=15 if p506r4>=7700 & p506r4<8400
	replace sector1r4=16 if p506r4>=8400 & p506r4<8500
	replace sector1r4=17 if p506r4>=8500 & p506r4<8600
	replace sector1r4=18 if p506r4>=8600 & p506r4<9000
	replace sector1r4=19 if p506r4>=9000 & p506r4<9400
	replace sector1r4=20 if p506r4>=9400 & p506r4<9700
	replace sector1r4=21 if p506r4>=9700 & p506r4<9900
	replace sector1r4=22 if p506r4==9900 

	#delimit ;
	lab def sector1r4 
	1 "Agricultura, ganadería y silvicultura" 
	2 "Pesca y acuicultura" 
	3 "Explotación de minas y canteras" 
	4 "Industrias manufactureras" 
	5 "Suministro de electricidad, gas, vapor y aire acondicionado"
	6 "Suministro de agua; evacuación de aguas residuales, gestión de desechos y descontaminación" 
	7 "Construcción" 
	8 "Comercio al por mayor y al por menor; reparación de vehículos automotores y motocicletas" 
	9 "Transporte y almacenamiento" 
	10 "Actividades de alojamiento y de servicio de comidas" 
	11 "Información y comunicaciones" 
	12 "Actividades financieras y de seguros" 
	13 "Actividades inmobiliarias"
	14 "Actividades profesionales, científicas y técnicas" 
	15 "Actividades de servicios administrativos y de apoyo" 
	16 "Administración pública y defensa; planes de seguridad social de afiliación obligatoria " 
	17 "Enseñanza" 
	18 "Actividades de atención de la salud humana y de asistencia social" 
	19 "Actividades artísticas, de entretenimiento y recreativas" 
	20 "Otras actividades de servicios"
	21 "Actividades de los hogares como empleadores; actividades no diferenciadas de los hogares como productores de bienes y servicios para uso propio"
	22 "Actividades de organizaciones y órganos extraterritoriales"; 
	#delimit cr
	lab val sector1r4 sector1r4

end
**************************************************************************
************************ENAHO*********************************************
**************************************************************************

use "$enahoData/enaho_2019-Modulo05/enaho01a-2019-500.dta", clear


append using "$enahoData/enaho_2020_tri1-Modulo05/enaho01a_2020_500.dta"
append using  "$enahoData/enaho_2020_tri2-Modulo05/enaho01a_2020_500.dta"
append using  "$enahoData/enaho_2020_tri3-Modulo05/enaho01a_2020_500.dta", generate(app_tri3)



addSectorLabels

replace aÑo = "2020" if app_tri3 == 1


gen realYear = aÑo
destring realYear, replace
destring mes, replace
gen _date = mdy(mes,1,realYear)
gen quarter = quarter(_date)
gen qdate = qofd(_date)
 

replace fac500 = fac500a if fac500 == . 

// qui: su fac500 if realYear == 2019 & quarter == 1, detail
// di "`r(sum)'"
// qui: su fac500 if realYear == 2019 & quarter == 2, detail
// di "`r(sum)'"
// qui: su fac500 if realYear == 2019 & quarter == 3, detail
// di "`r(sum)'"
// qui: su fac500 if realYear == 2019 & quarter == 4, detail
// di "`r(sum)'"
// qui: su fac500 if realYear == 2020 & quarter == 1, detail
// di "`r(sum)'"
// qui: su fac500 if realYear == 2020 & quarter == 2, detail
// di "`r(sum)'"
// qui: su fac500 if realYear == 2020 & quarter == 3, detail
// di "`r(sum)'"


/*----------------------------------------
-- Need to reweight to make sure all quarters hacve same weight 
----------------------------------------*/


bys qdate: egen double wgtsum_quarter = total( fac500) if realYear == 2019
bys realYear: egen double wgtsum_2019 = total( fac500) if realYear == 2019

*multiplication factor - scales up a quarter's weight so it returns the 2019 annual total. 
gen double scale_factor = wgtsum_2019 / wgtsum_quarter


gen  double  test_weight = fac500 * scale_factor





// the main problem is 2019, I need to scale up each quarter so we get the 2019 total.
qui: su test_we if realYear == 2019 & quarter == 1, detail
di "`r(sum)'"
qui: su test_we if realYear == 2019 & quarter == 2, detail
di "`r(sum)'"
qui: su test_we if realYear == 2019 & quarter == 3, detail
di "`r(sum)'"
qui: su test_we if realYear == 2019 & quarter == 4, detail
di "`r(sum)'"




replace fac500 = test_we if realYear == 2019

***** 
keep ubigeo ocu500 fac500 sector1r4 p507 p523 p524a1 p530a p538e1 conglome estrato quarter realYear qdate
//keep employed
keep if ocu500==1


***Nivel dpto
gen dpto=substr(ubigeo,1,2)
order dpto




/*----------------------------------------
-- Collapsing with weights
----------------------------------------*/


collapse (count) N = fac500 (sum) weighted_N_national =  fac500 ,by(sector1r4 qdate realYear quarter)



save "$data_folder/sector_enaho_2020_2019_quarterly_right_weights.dta", replace





